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Abstract 

We study the Hawking radiation for the geometry of an evaporating 1+1 dimensional 
black hole. We compute Bogoliubov coefficients and the stress tensor. We use a recent 
result of Srednicki to estimate the entropy of entanglement produced in the evaporation 
process, for the 1+1 dimensional hole and for the 3+1 dimensional hole. It is found 
that the one space dimensional result of Srednicki is the pertinent one to use, in both 
cases. 
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1. Introduction 



Several interesting phenomena are related to the discovery of Hawking radiation 
It is intriguing that black holes seem to obey laws of thermodynamics 0. The 
information contained in the matter which made up the black hole is lost into the 
singularity. Hawking radiation appears in the evaporation of the hole, but the outgoing 
modes are not in a pure state; instead they are mixed with modes of the field that fall 
into the singularity. The precise significance of black hole thermodynamics, and its 
relation to the ordinary ideas of thermodynamics and information theory, are matters 
of debate. 

Recently the discovery of 1+1 dimensional models for black holes [§, ^ has led to a 
more accurate understanding of the semiclassical features of black hole geometry and 
Hawking radiation. In particular the model of |Q (the RST model) can be exactly 
solved to yield the semiclassical geometry of a black hole formed by a shock wave of 
infalling matter, and evaporating by massless scalars to an endpoint (the 'thunderpop'). 
It may even be possible to obtain a complete quantum gravity plus matter description 
of the black hole evaporation process in 1+1 dimensions 

In this paper we study some features of the semiclassical geometry and Hawking 
radiation in semiclassical models. For the RST model of the evaporating hole we com- 
pute the Bogoliubov coefficients. We perform a point splitting calculation to compute 
the stress tensor at T+. We also compute the stress tensor in the evaporating region 
using the anomalous trace of the matter stress tensor. The RST model is solved also 
for the case where the hole is formed by one shock wave, evaporation occurs for a time, 
and then a second shock wave increases the mass of the hole again. (This geometry 
is used for clarifying some aspects of the entropy produced by the hole, as discussed 
below.) 

We then study the 'entropy of entanglement' of the Hawking radiation, by two 
methods. We can compute the density matrix obtained by tracing the field modes 
inside the horizon. This was the approach taken by Hawking and also carried out in 
for the 1+1 dimensional case, and it yields a density matrix that is close to thermal after 
the initial stage of the collapse and formation of the hole. We consider the case where 
the hole has a finite lifetime (due to the evaporation) and thereby estimate the entropy 
of the entire radiation produced. We then compare this result to that obtained by using 
a calculation of Srednicki 0. Srednicki considered a scalar field in fiat Minkowski space 
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in the vacuum state, and 'traced out' the degrees of freedom inside a ball of radius R. 
The entropy of the reduced density matrix is the 'entropy of entanglement' between 
the region inside the ball and its complement. The entropy S depends on R and also 
on the ultraviolet cutoff, which gives the 'sharpness of separation' between the regions. 
In the one space-dimension case, the infrared cutoff also appears. We find that both 
for the 1+1 dimensional black hole and for the 3+1 dimensional hole the one space 
dimension result of Srednicki is the pertinent one to use, and the leading dependence 
of the entropy on the black hole mass is reproduced. 

The plan of this paper is as follows. In section 2 the RST model is reviewed. In 
section 3 the Bogoliubov coefficients for a scalar field in the evaporating black hole 
background are computed. Section 4 contains the calculation of stress-tensor. Section 
5 studies the two shock wave solution. We discuss the entropy of the Hawking radiation 
for 1+1 dimensional black holes in section 6, and for 3+1 dimensional black holes in 
section 7. Finally, a discussion is presented in section 8. 

2. The RST model 

The model of Russo, Susskind and Thorlacius (RST) [|] is a modified version of 
the model of two dimensional dilaton gravity coupled to quantum matter introduced 
in . The key idea of RST was to introduce an additional counterterm which restores 
a global symmetry originally present in the classical dilaton gravity + matter action. 
This allowed them to solve the theory analytically in the large limit. The properties 
of the RST model have been extensively discussed in the literature 0, ||, |^ so we will 
just mention the facts we will need for later use. 

The semiclassical effective action of RST can be written as follows 

^ r N 1 ^ 

- (fx^-g{x) j (fx' ^ -g{x')R{x)G{x,x')R{x') , 

where R is the 1+1 dimensional scalar curvature, is the dilaton field and /j, i = 
1 . . . ,N are N massless conformal scalar fields. G{x, x') is the Green's function for the 
d'Alembertian in curved space. The constant A plays the role of Planck mass. 

The analysis of the semiclassical equations of motion that follow from the action 
(1) can be simplified by the following two steps. First, one can write the metric in the 



2 



conformal gauge, given by 

where = x° ± x^. Secondly, one can make a field redefinition and introduce the 
fields 

n = «:-e-^ + ^ + lln^ (2) 
X = + p - ^ - ^ ln(4«:) , 

where k = The coordinates can be fixed so that = so the dilaton field 
differs from the conformal factor of the metric p by a constant. The matter fields 
are assumed to reflect from the strong coupling boundary VL = Qcr = 1 0- The 
semiclassical equations can now be reduced to 

d+d^f, = (3) 

Here are the components for outgoing and ingoing matter energy of the stress- 
tensor, which is defined as follows. Since the classical matter action is written as 

N 

I (i^.T^/ — O 

An 



Sf = -^ [d'xV^Y.i^0 , (4) 

i=l 



the stress tensor is 

f ^ -2 SSf 

The components representing outgoing and ingoing matter are normalized as 

271 . 

in the conformal gauge. 

The functions t±{x^) are fixed by boundary conditions. We assume that the in- 
coming matter energy flux at I~ vanishes sufficiently rapidly at early and late times. 
Then, t+(x+) = l/4(x'^)^. The solution for the field Q can now be found to be 

n = -X^x+(x- + 4^P+(x+)) + ^M{x+) - \ ln(-4A'a;+x-) , (5) 
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where 

M(x+) = X r duu Ti+{u) (6) 
Jo 

p+{x+) = r du tiau) . 

Jo 

Consider now an incoming matter shock wave that carries energy M. The stress tensor 
is then given by 

tUx^) = ^6{x^ ~ xt) , (7) 

which is the only non-vanishing component. We substitute this in the equations (6) 
above. Following we set Xxq = 1. 

After solving the RST equations one finds the following results. The Penrose dia- 
gram for the evaporating black hole spacetime is given in Fig. 1. 

The spacetime is seen to be divided into three regions. The 'lowest' region I, before 
the incoming shockwave at Xq, is the linear dilaton vacuum, bounded from the left by 
the timelike strong coupling boundary Q = Qcr- The incoming shockwave forms the 
black hole by causing the boundary Q = Q^r to become spacelike. It can be shown that 
the scalar curvature R diverges at the spacelike portion of the boundary. The apparent 
horizon also forms after the incoming shockwave. Following the apparent horizon 
is defined by the condition 9+0 = 0. After the black hole forms, it starts to evaporate, 
and the apparent horizon shrinks until it meets the singularity at the endpoint of 
evaporation. At the endpoint of evaporation a short (delta function) burst of negative 
energy is seen to emerge from the black hole. This is called the 'thunderpop' [^, and it 
travels along the null line x~ = x~ to future null infinity. ( IS cL light-cone coordinate 
of the endpoint to be specified later.) Thus the region II bounded by the thunderpop 
and the incoming shockwave is the curved region of the evaporating black hole. After 
the thunderpop, the spacetime becomes again fiat and the boundary Q = fl^r becomes 
timelike. The corresponding region III is a linear dilaton vacuum. 

In the linear dilaton vacuum region I the metric is 

ds'^ = {X'^ x^ x^ dx'^ dx~ . 
We can write it as ds"^ = —dy^dy" using coordinates 

y~ = -^ln(-Aa;-) (8) 
y+ = ^ ln(Ax+) - y^ . 
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The shift is introduced to set the origin of the coordinate to the point A where 
the reflected ray of the event horizon (see Fig. 1.) meets X~ . 

The event horizon, the singularity and the apparent horizon meet at point E, the 
end point of evaporation. In our conventions, its coordinates are 

fe<) = (-^(l-e-^)-', jj5(e^-l)). (9) 

We can now specify what the shift is. In region I the reflecting boundary Q = Qcr- — 
\ is the line 

+ -1,1 + 

y '^Y^^^y^- 

Reflecting the line x~ = xj (off the boundary Q = Qcr) back to X~, we find that the 
point A has y^ = —j ln(— Aa:~) + ^ In \—yQ. Setting y~^ = yields yQ = —j ln(— 4Aa;~). 

The region II is the curved region of the evaporating black hole. It can be joined 
continuously but not smoothly with region III along the line = , with region III 
being again a linear dilaton vacuum, but with the coordinate x~ shifted. From the 
joining conditions the metric in III can be found to be 

ds^ = (X''x+(x- + ^))-^dx+dx- . 

If one makes the coordinate transformation 

(7+ = \\n{\x+) (10) 
A 

this metric becomes ds^ = —da^da~. The coordinate a~ has been normalized so that 
the thunderpop is at = 0. The reflecting boundary Q = Qcr in region III is at 
a+ = a- + \HXxt). 

The metric in II becomes asymptotically flat near X+. We can extend the coordi- 
nates (T^ from region III into region II in the neighbourhood of X+. Then the metric 
in region II also has the asymptotic form ds^ — > —da~^da~ near X+. Thus, are the 
physical coordinates near X+. 

Finally, we identify some points of interest in the Penrose diagram. The point A 
where the reflection of the null line x^ = xj meets X^ we already set to be at y^ = 0. 
The point B is the projection of the end point E along a null ray to past null infinity. 
We find it to be located at y^ — ^0-- The point C is defined by projecting the point 
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where the apparent horizon meets the incoming shockwave along a null ray to future 
null infinity. It is at cr^ = ln(^(exp(^) - 1)). Since M » nX, to a good 
accuracy = — The absolute value of is thus the total (retarded) time of 
evaporation of the black hole. 

3. Bogoliubov transformations 

In this section we calculate the Bogoliubov coefficients for the relation between the 
natural mode functions in the 'in' region close to X~ and the 'out' region close to X~^. In 
||l[] the Bogoliubov coefficients were estimated for modes travelling close to the horizon 
that forms in a spherical collapse of a star. In [^] the Bogoliubov coefficients were 
computed for the 1+1 dimensional eternal black hole geometry, i.e. without taking 
into account the backreaction on the metric due to the Hawking evaporation. Our 
notations follow those in |0 where we refer to for all introductory steps. 

We choose the following form for modes at X^, X+ respectively 

= -^e-'^y^ (tn) (11) 
\/2uj 



^ (out) 



'2u 

where the normalization factor is in agreement with the form (^ of the matter action. 

We define the Bogoliubov coefficients a^^^)' and (3u]u)' for the relation between the 
'in' and 'out' modes as follows 

/■oo 

V^= duj' [a^^'Mo;' + Pcouj'Ut'] ■ (12) 

Jo 

To compute the Bogoliubov coefficients, we pull the 'out' mode back to X^. First 
we divide the mode into two parts at the point cr^ = 0. The 'upper' piece v^j = 



/2w 



e ^^'^ 6 {a ) reflects from the timelike boundary in region III without experiencing 



any blueshift. At X it becomes 



1 e-M.+ -4)^(^+_^+) 



^2uj 

where was defined in the previous section. 

However, the 'lower' piece = g^ts distorted. It first experiences 



a blueshift when pulled back to region I. This is done by replacing the coordinate a 
with the coordinate y~ using the relation 

^ ~ ~X Xx~ + tiM/Xk ^ ■ 
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Then we reflect the mode from the timehke boundary in region I back to X , by 
replacing y~ with 

= + yt -\'^'^\ = -\ in(-Ax;) . 

The final relation between the coordinates a~ and can then be written as 

= -lln[AA(e-^^^-C)] , 



where 



AA = = e^-^/^'^ (14) 

Ax- + ttM/Xk ^ ' 

_ -47rM/AK 

Axr 



As an important aside, notice that (|13]) implies a relation between a small distance 
da^ at X+ centered at a~ and the corresponding small distance dy^ at which results 
from mapping the former distance along null rays which reflect from the boundary back 
to the past null infinity. This relation can be found to be 

dy+ = {1 + (e^-M/.x _ ^ya- _ (^5) 

The result (^) tells us that a distance da^ at X+ (before the endpoint of evaporation) 
maps to a much smaller distance dy~^ at 2~ , with the 'squeeze factor' becoming ex- 
ponentially larger as the black hole evaporates. This observation will turn out to be 
crucial in applying the Srednicki calculation for black holes, as will be discussed in the 
section 6. 

Returning back to the behaviour of the modes, ([T3|) tells us that the 'lower' part of 
becomes 

= ^exp{^ln[AA(e-^^^ -C)]}^(-1/+) (16) 



V2u A 
when pulled back to X^ . 

It is now straightforward to proceed to find the Bogoliubov coefficients. The result 

is 



1 luj' r /-o 



{ f dy+ exp{^ ln[AA(e-^^^ - C)] + lu'y+j (17) 

J — CO A 
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^e'^^^B I dz+expi-ifcu-cu')'^^}} 

^0 

/3u,u,' = ^XHif dy^exp{'-^ln[XA{e-'y^ -C)]-zu;'y^} 

ITI \ UJ J-oo A 

■ I + r°° 

+ e-"^ / dz+ exp{-i(a; + u')z+}} . 
Jo 

These results resemble the ones of Giddings and Nelson (GN), with the following three 
differences: 1) there are additional terms resulting from the 'upper' part of the mode 
Vi^, i.e., the part after the endpoint of evaporation, 2) AA is different, and 3) we have 
C = while GN had C = 1. (Of course, C ^ 1 since M/kA » 1. However, 

the small difference turns out to be significant if one tries to investigate the behaviour 
of the modes very near the endpoint.) 

The integrals in the above expressions can be evaluated further. Substituting first 
X — e^y^ and then t — Cx, we get 

2tt uj Jo 

^^^'vt — z! — 1 . 

iij — Lo' — ie 

The remaining integral can be identified as an Incomplete Beta function. The coefficient 
/5a;a;' IS Computed similarily. The final expressions are 

^ ''^')VM(AAr/"C^('^+'^')/^5c(-^ + ^ + 6,z + ^) (19) 



27rA Lo AAA 

13..' = ^ i-y^' {(AA)-A c^i^-^'yx Bci-'-^ _ !^ + e, , + !^) 

ZnA UJ AAA 

-iX e-"^'yB + oj' - ie')-^} . 

Note that in the semiclassical approximation the thunderpop is a delta function at X"*", 
thus there is a part of the field modes that is not captured by the modes e^*'^°" at X"*" 
for any finite range of uj. Thus the Bogoliubov coefficients computed below need to be 
supplemented with an infinite frequency component to completely describe the field at 
X+. 

Let us now turn to the issue of examining the nature of the radiation by studying 
the approximate behaviour of the Bogoliubov coefficients. We expect outgoing thermal 
radiation at constant temperature ^ to be seen in the region cr" G (~~^)0) of 
X+, except perhaps at the beginning and very end of the evaporation process. The 
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coordinate corresponding to this region is very small, so we can approximate 

ln[AA(e-^^'' - C)] ^ ln[l - e'"^'/"^Ay+] . (20) 

If we are not too close to the endpoint, the term 1 in ( pO]) is negligible. For the essence 
of Hawking radiation, we have the frequency ~ A at I'^ and very high frequencies 
ijj' ~ ;^e47r(A/+j/+)/KA --£- such 10 ^uj' wc cau ignore the second integrals in the 

expression ([171) for a^^^/ , fi^^^l . We therefore get 

^ ^^{f rfy+expl^lnhe^-^^A'^Ay+l+^cV}} (21) 

I'K \ UJ J~oo A 

Zn \ uj j-oo A 
This means that we get the thermal relation (see 0) 

a^.' ~ -e""/"/3..' (22) 

and we see that the outgoing radiation is thermal with constant temperature Th = 
in the region a~ e (-^, 0) at J+. 



4. The stress tensor 

Since we have worked out the relations between the 'in' and 'out' modes, we can 
easily calculate the VEV of the stress-energy at 2'^, i.e. {0,in \ T^y{a^) \ 0,inY^'^ by 

using the point-splitting method. The non- vanishing component of {T^j^uY'^"' is the 

component. Since we are computing the 'in' VEV at X+, we need first the form of the 
'in' mode at X'^. For cr^ < it is 

= -|=exp{^ln[^e-^- + Cm-a~) . (23) 



Again, for (T~ > (after the endpoint) there is no redshift and it is then trivial to see 
that the stress-tensor vanishes in this region. We concentrate only in the region before 
the endpoint. Using the point-splitting method, we first calculate 

1 (9 (9 

(T__(a-))= hm ----- dcu u^{a^)uUa^) , (24) 



*cr- 2tt dai d(72 "'0 

where the LHS means the VEV with respect to the 'in'-vacuum. (We use this notation 
in the following.) For a" < the step functions are just constant and can be ignored. 
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Taking the partial derivatives and making a series expansion in (af — aol) yields then 
a term which diverges quadratically in the limit, and a finite term. The divergent term 
must be subtracted (it is just the usual vacuum divergence) and the finite term gives 
the renormalized value for the stress-energy: 

1 

(T_(a-))^'^" = —[1 - n +fe4-^/M-l)e^-)2] • (25) 



This result is for one conformal scalar field, for fields it must be multiplied by 
N. Note that in the region E approximate 

(^__(a-)r«^^ = ^(^^)^ (26) 

which is the correct value for outgoing thermal radiation at temperature Tjj = 

(Comment: if one would use the approximate behaviour of the mode very near the 

end point, one would find that (T Y'^'^ ^ at roughly Planck distance from the 

end point. One should not, however, trust such a local treatment when dealing with 
modes.) 

The above formula gave (T ) only at X~^. However, in 1+1 dimensions it is easy 

to calculate {T^^Y^"^ everywhere by using the trace anomaly, integrating the covariant 
conservation equation and applying the boundary conditions (see for details). Us- 
ing this route, we end up with the following results for the stress tensor everywhere in 
region II: 

~ 487r(a;-)2^e-2p-l/4^^^ e-2p-l/4 ^ ^ 

^ N 1 e-^^ + 1/4 A^x+(x-+7rM/^A^) + l/4 , 
^ 487r(a;+)2^e-2p-l/4^^^ e-2p-l/4 ' ^ 

yen ^ N \' 1 6'^^ [A^x^ (x' + ^) + 1/4] [A^a:^^- + 1/4] 

^ 247r e-2p- 1/4^ (e-2p- 1/4)2 A2x+x- ^' 

where p is given implicitely via the equation 

' + S = -A'x+ X- + -—) - - In -A^'x+x- + ^ . 28 
/ kA^ 4 k\ 

In the above the stress tensor is given in the 'Kruskal' coordinates x^ since one is 
generally interested in its behaviour in both sides of the apparent horizon. 

5. Two incoming shock waves 
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In this section we present the solution of the RST equations in the case of two 
incoming shock waves. We use this geometry later in section 6 where we discuss the 
entropy of the black hole. We want to first form a black hole of mass Mq, which then 
starts to evaporate. Then the second shock wave at a later time carries additional 
energy to the black hole. For instance, it could restore the black hole back to its 
original mass. The question then is: How much is the entanglement entropy after the 
evaporation process? If the second shock wave carries just enough energy to restore 
the black hole back to its original mass, but not more, is this entanglement entropy 
related to Bekenstein entropy? We discuss this question in section 6, here we will just 
derive the results that we will use. 

Recalling how we defined the incoming shock wave, it is clear that for two incoming 
shock waves we should replace (|^) with 

tUx^) = ^6{x^ - 4) + - xt) . (29) 

The first term on the RHS is the first shock wave at Xq (we set again Xxq = 1) carrying 
energy Mq and the second term is a later shock wave at xf carrying energy Mi. We 
do not specify xf, Mi here, but xf should be chosen so that the black hole formed by 
the first shock wave has not yet evaporated away when the second shock wave reaches 
it. On could think of Mi being the energy needed to restore the black hole back to its 
original size, but for now that is not essential. 

The spacetime curve for the apparent horizon can be found by solving the equation 
= (this follows from setting 9+0 = for the dilaton field; which is the RST 
definition of the apparent horizon). The solution is 

^ ~4A(x- + ^ + Jft^^(x+-x+)) ' 

for x^ > Xq. The boundary of the spacetime is the critical line Q = Qcr = \, which 
defines another curve x~^ = x~^{x~). The (final) endpoint of evaporation is where the 
above two curves meet. We find it to be located at 



X 



TT Ml 



(Mo + T^)(l - e-^"(*^"+^^i)/"^)-i (31) 



kA^ Xxi 



X 



_^g47r(Mo+Mi)/KA _ 
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We can now join the region II with the linear dilaton vacuum in region III with an 
appropriate shift of the coordinate x~ . For the metric in region III, we find 

as 



This in turn tells us how to define the coordinates which are the 'physical' coordi- 
nates near . We define 

(j+ = ^ln(Ax+) (32) 
A 



a- = --ln[ 



A '\x- + f^{MQ + Mi/\x\, 
On the other hand, we still have the 'physical' coordinates in region I 

y+ = \H\x^)-yt (33) 
y- = -^ln(-A2;~) , 

where yQ = — ^ln(— 4Ax7), with x~ given now by (|3T|) . The reflecting boundary in 
region I is the line y~^ = y~ + jln ^ — yQ . 

It is straightforward to see that the relation between the coordinates , y^ now 
becomes 

a- = -\ ln[AA'(e^^^^ - C')] , (34) 
A 

where 

AA' = e^"(*^°+^^^)/^" (35) 

(J' = I — g-47r(Afo+Mi)/AK 

As we noticed in section 3, this implies a relation between a distance da^ at X+, 
centered at a~ and the corresponding small distance dy^ at X~ . In the two shock 
wave case, the relation is 

dy+ = {l + (e4-(A^o+Mi)/«A _ ^yXa' yl^^- ^gg^ 

Again, a distance da~ in X"*" (before the endpoint of evaporation) maps to a much 
smaller distance dy~^ in X~ and this 'squeeze factor' becomes exponentially large as 
the black hole evaporates. In the present case the 'squeeze factor' reaches the value 
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g47r(Mo+Mi)/KA ^hich exceeds the value ~ q'^-^Mo/>^>^ that would be obtained in the 



absence of the second shock wave. The result (|36D will be used in the next section. 



6. Entropy for 1+1 dimensional black holes 

As argued by Hawking, the process of pair creation by the gravitational field of the 
black hole creates a state which is 'mixed' between the regions exterior and interior to 
the horizon. If we compute the reduced density matrix that describes the field in the 
exterior region, then the entropy computed for this density matrix gives the 'entropy 



of entanglement' [|lT| between the interior and exterior regions of the hole. If we do not 
take into account the backreaction from the created radiation then an infinite number 
of particle pairs are produced and the entropy of entanglement will be found to be 
infinite also. But in the simple RST model we can estimate the produced entropy for 
the semiclassical geometry that includes backreaction. We shall do this in two ways: 

1. We directly compute 

S = -Tr{p\Yip} (37) 

where p is the reduced density matrix describing the field in the exterior region. 
Here we use the fact that p is to a good approximation a thermal density matrix. 

2. We use the result of Srednicki described in the introduction. Thus a complete 
spacelike hypersurface is considered, and the part corresponding to the interior 
of the horizon is assumed to be, effectively, the traced over region considered in 
this approach. 

1. The essential idea is to define reasonably localised regions on such that the 
density matrix can be described as approximately thermal in those regions. Hawking 
presented this analysis for the 3+1 dimensional black hole |]l|; it was worked out ex- 
plicitly for the 1+1 case in (without backreaction). We follow the notations of |^ 
in the following. Define a complete set of orthonormal wavepackets on X"*" : 

1 in r^j+'^y n ■ I 

v^n = / du e^'^*""/^ , (38) 



where j = 0,1,2,... and n are integers. The wavepacket vjn is peaked about X"*" 
coordinate a' = 27m/ e, has a spatial width ~ and a frequency uj ^ je. In this 
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basis the reduced density matrix obtained by tracing out the field inside the black hole 
is 

p = N-'J2 e-^^-"^""1KJ)({n,„}| . (39) 

{rijn} 

This has the form of a thermal density matrix. We can write 

P = IlP.n (40) 



with 



Pjn = Nj^n r{njn)\njn){njn\ (41) 

njn=0 

r{njn) = e"^"^"'^^- . (42) 



We get 



where 



S = j:S,n (43) 



Sjn = -Tr(^jn){Pjn^^Pjn} = - r{njn)liar{njn) ■ (44) 

njn=0 

A brief computation gives 

Sjn = pUjie'^'^^ - 1)-^ - ln(l-e-^'^0- (45) 

Thus 

oo oo 

= E^.n = Yil^^A^^^' - - ln(l-e-^-0} (46) 

j=0 j=0 



dj{pje{e^^' -ly^ - ln(l -e"^^'^)} 



/o ' " ' ' ' " 3/5e 

The entropy from evaporating fields is 

S = NY.Sn. (47) 

n 

The separation between wavepackets ( pSj ) is Acr^ = 27r/e. Thus in time T the number 
n ranges from n = 1 to n = tT /2t{. From (^, we get 
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The total evaporation time is T = ^f^- This gives the estimate of the total entropy 
created in the Hawking radiation: 

47iM , , 

5total = . (49) 

Note that this is twice the Bekenstein entropy, which for the 1+1 dimensional hole 
is S'sek = The result ( ^91) is the entropy of a thermal distribution of bosons at 



temperature /3 ^ and with energy E = M, which in one dimension is given by S* = 2i3E. 
For a discussion of how the entropy of radiation at X+ relates to the Bekenstein entropy, 
see 0. 

2. We now investigate the application of Srednicki's result to the black hole. We 
split the discussion into 3 parts: 

• (i) We recall the result of |^, and discuss the issue of infrared divergence. 

• (ii) We discuss how this result may be applied to the black hole, after making 
some plausible arguments for the required modifications. 

• (iii) We carry out the calculations for the entropy. 

(i) The computation of Srednicki for the one space dimension case may be described 
as follows.[] We consider a free scalar field on a 1-dimensional lattice, with lattice 
spacing a. Let this field be in the vacuum state. We select a region of length R of 
this lattice and trace over the field degrees of freedom outside this region. This gives 
a reduced density matrix p, from which we compute S = —Tr{plnp} which is the 
entropy of entanglement of the selected region with the remainder of the lattice. It is 
immaterial whether we trace over the interior or exterior of the selected region; since 
the field was in a pure state overall the entanglement entropy is the same in both cases. 
This entropy is given by 0] 

S = Ki\n{R/a) + K2\n{fiR) , (50) 

for one scalar field. (For species the result must be multiplied by A^.) Here n is an 
infrared cutoff. 

The infrared term is very sensitive to boundary conditions. As an example consider 
taking a periodic lattice, and let the scalar field be periodic as well. If the field is 



Similar issues were studied in [[L3l 
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massless, then the zero mode of the field varies over an infinite range, in the vacuum 
state. If we trace over the degrees of freedom in a subregion of the lattice, then the 
mean value of the field inside is correlated to the mean value outside, but this mean 
value can take on an infinite range of values. The entanglement entropy will thus be 
infinite. 

If we take antiperiodic boundary conditions for the scalar field then the zero mode 



does not exist. With this choice (Ew becomes 14 



S = l-HR/a) + ]:HI/R) (51) 

where / is the infrared cutoff coming from the finite size of the lattice. Another way 
to kill the zero mode of the scalar field is to have a vanishing boundary condition for 
the field at say x = 0. Let the traced over region extend from x = Xi to x = X2. 
Now modes with wavelength much greater than X2 effectively vanish over the interval 
(xi,X2), so they do not serve to entangle this region with the remainder of the line 
X > 0. Thus the entanglement entropy will be finite without the need for an explicit 
infrared regulator. A third way of dealing with the zero mode is to simply assume that 
the field has a small mass /i. Then no other infrared regulation should be needed and 
the results should not be sensitive to choice of boundary conditions. 

(ii) For the black hole, we start by considering a complete spacelike hypersurface 
through the evaporating geometry, described as follows. Starting at spatial infinity 
(cr^ = — oo), we move near X+ to a point with cr^ = < 0. (The black hole vanishes 
at = 0.) Then we smoothly bend this hypersurface so that it enters the horizon and 
reaches the timelike segment of the line Q = Qcr- (This timelike segment occurs before 
this critical line becomes the singularity.) Thus the hypersurface is kept spacelike 
all through, and avoids the singularity by passing below it. An observer collecting 
radiation far from the black hole will see the part of this hypersurface that lies along 
T"*", and we would wish to trace over the part that cannot be observed from outside 
the black hole. This would give the reduced density matrix, and thereby the entropy 
of entanglement between the field inside and outside the hole. Let us be more precise 
about what we take as the 'observed' part. Suppose that the observer at T"*" carries 
an instrument with her which she uses to collect the outgoing radiation. At first, near 
a~ = — oo nothing comes out from the black hole. The observer has to wait for quite 
a while before the black hole starts to radiate. At some point the radiation starts 
and becomes approximately thermal. We have discussed earlier that this happens 
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roughly at = Correspondingly, around this point the observer turns on her 

instrument. The observer then collects radiation up to some retarded time , when 
she turns her instrument off again. Thus the part of the hypersurface between and 
(jf corresponds to observations, and the rest of it will be traced over. Notice that the 
observer can neither start nor stop collecting radiation at an instant, but there will be 
a short time scale da~ which she needs to turn on or shut off her instrument in a proper 
way. The time interval needed should be sufficiently rapid to give a good accuracy for 
specifying the turn-on point and the shut-off point , on the other hand it should 
not be so small that it creates disturbances in the matter field that generate radiation 
comparable to the Hawking radiation collected. It is reasonable to assume that the 
time intervals d(T~{aQ ), da~{ai ) needed should be given by the typical wavelength in 
the outgoing thermal radiation. Thus we can assume that da~ ~ (3h ~ \ (our result 
for the entropy of the hole will turn out not to depend on this choice). 

We split the contribution of different modes to the entropy, as follows: a) the modes 
of wavelength cj"^ >> Mj}?: these modes are effectively constant over the observed 
region ((7(J",(7f), so they may be taken as a contribution to the 'zero mode'; b) the 
leftmoving modes with cu"^ ~ M/A^; c) the rightmoving modes with o;"^ ~ M/\^. 

Modes of type a) will give a divergence in the entanglement entropy even in flat 
Minkowski space (without black hole). This happens when the mass ji of the field is 
taken to zero, or (if the field is massless) as the observed part of the hypersurface is 
taken further and further away from the line Q = Q^cr where the field vanishes. Since 
we are interested in entropy of entanglement of the rightmoving Hawking radiation 
(which occurs over a distance M/A^), we subtract this (diverging) contribution arising 
from the large wavelength modes » M/\^. We also ignore the leftmoving modes 
b) , as they do not contribute to the Hawking radiation. The rightmoving modes c) are 
of interest to us, but after particle pairs have been created, the state of the field is not 
the vacuum state for the geometry on the spacelike hypersurface under consideration. 
Srednicki's result, on the other hand, apphes to a vacuum state for the field. The 
essential idea is to follow the rightmoving [i.e. outgoing) modes of the field back from 
T+, through the collapsing matter to the line Q = Qcr- where they reflect to left moving 
modes that can be followed back to X~ . Here these left movers are in the vacuum state, 
so that we may apply a Srednicki type approach to estimate the entanglement entropy 
in these modes. 

As we follow the radiation modes back to X~ in the manner indicated above, we 
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observe the following. The region G {a^ ^a^) in X+ corresponds to a region G 
iy^ {o-q) , {ai)) in X", where the relation between a~ ,i = 0, 1 and y^{al),i = 0, 1 is 
given by the formula (|13D in section 3. Thus, the latter region is the region of starting 
points for ingoing rays which will experience redshift and give arise to the collected 
radiation in the 'observed' region of Also, this region is separated from the rest 
of X~ by 'cuts' of size dy^{aQ) and (i?/+(o"f ), the size of which follows from the size of 
the corresponding cuts near X"*" by the relation ([ISD given in section 3. Now we can 
apply the result of 0. We disentangle the finite region y^ G (?/'^((7o ), ?/+((Tf )) from 
the rest of X~ by cuts of size = (i?/+(cro ), cti = dy^^a^). Ignoring 'zero modes' and 
the leftmovers, we expect to create an entropy for each scalar field (see discussion in 
section 8 below) 

(52) 



R R 
S = K3 ln(— ) + /t3ln(— ) 



where R = ?/^(o"r) ~ V^i^o) ^^^^ f^s = i^i- (A factor of | because we are considering 
only rightmovers and another factor of | because the contribution to S is separated 
over the two 'cuts'.) 

"™ as 



fiii) Let us rewrite 



2/t3ln(i?A) + ln( 



1 



^3 ln(' 



1 



We substitute aj = dy^ 



ai 



0, 1 given by the relation ([151) : 

{1 + (e4.A//KA _ 
-AttM/kX -Xa^ ^ 



l)e^''^}-^da~ 



(53) 



(54) 



A 



1 
A 



The former approximate formula is valid for ai G {—^0-, 0) and the latter just results 
from the fact that the redshift is negligible for the rays at earlier times. 



After substitution we get (ignoring the infrared cut-off term) 

S = fi:3ln(e^^^/''V'^i") + /€3ln(YA) + 2K3ln(i?A) 

A 

AtiM 

= ^3{^^ + (xr) + 2/t3 ln(i?A) . 



(55) 



The second term is negligible with respect to the first term, since it can be calculated 



that R 



The first term is the significant term. As we take closer and closer 
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to the endpoint a =0 (the observer collects more outgoing radiation), the first term 
approaches 

. (56) 

The above result for entropy must be multiplied by A^, the number of scalar fields. 
Let us deduce the value of K3 by comparing (^) with the entropy of entanglement 



estimated directly from the density matrix of the outgoing radiation (eq. (|49|) ). This 
gives 

5-^,if.3 = -. (57) 

If our assumptions regarding the separation of the 'zero mode', of right and left movers 
are correct {i.e., if ki = 4^3, ki as given in (^)), then ( p7D agrees with the calculation 
of Srednicki which gave ki = |. We discuss this issue further in section 8. 

Note that if we consider the evaporation right up to the endpoint, then the cutoff 
scale (over which the radiation measurement is switched off) must go to zero, and 
the entropy becomes infinite. But since we are using the semiclassical geometry, it 
is not justified to go below distance (or (A^A)~^, for large A^) in our analysis. 
With this restriction, the entropy from the cutoff scale of the endpoint is ignorable 
compared to the entropy of the hole, for holes that evaporate over classical time scales 
A-kM/kX^ » \-\ 

The significance of the Bekenstein entropy for a black hole is a matter of debate. 
One hypothesis is that the horizon behaves as a membrane with e^^"'^ states, so that 
there is an upper limit to the entanglement entropy of the matter outside the hole 
with the hole itself [|I^]. Thus if a sufficiently large amount of matter were thrown into 
the hole then a part of the information would have to leak back out through subtle 
correlations in the Hawking radiation [1TB . 



In a semiclassical treatment of the gravitational field, on the other hand, there 
seems to be no limit to the amount of information that can disappear into the black 
hole. Thus the entanglement entropy can also grow without bound when matter is 
repeatedly thrown into the hole and the black hole mass allowed to decrease back to 
M by evaporation.^ It is possible to verify in the simple 1+1 dimensional evaporating 
RST solution that the entanglement entropy can indeed exceed the Bekenstein value 
by an arbitrary amount. We illustrate this by taking a simple example: the RST model 
with two incoming shock waves, discussed in the previous section. 



^We are grateful to S. Trivedi for pointing this out to us. 
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We again apply the Srednicki approach to estimate entropy. Consider an observer 
at X+ collecting radiation, who switches on the measuring device at aQ and switches it 
off at (jf with the corresponding switch-on-off intervals da^{aj^) as before. The only 
difference in the two shockwave case is that now the relation between the da~^s at 
and the corresponding dy'^^s at X~ is different. The relation was given by the formula 
(EB) in section 4. Now we need to substitute 



ai 



ao 



= {l + (e4-(^^«+*^i)/«^-l)e^"i"}-M(T- (58) 

^ -47r{Mo+Mi)/KA -Ao-j7 }_ 

A ' 

1 



A 

into the equation ( p3D above. Also the distance R will be different, but it is still of 
the order of ^ and the i?-dependent term can thus be ignored. All this gives for the 
entropy (for fields) 

S ^ A^/«3ln(e^"(^o+^i)/'^^e^"i") + ArK3ln(iA) + ... (59) 

A 

47r(Mo + Mi) _ 

^ Nk3{ + ai +... , 

kA 

where + . . . represents the ignored contributions from the subleading i?-dependent 
term and the infrared cutoff term. As ^0, the entropy becomes 

S^N.,'-^^Ml±^. (60) 

KA 

Substituting = 1/12 gives then 

^^MA^+MO^ (61) 
A 

Note that the Bekenstein entropy, on the other hand, need not exceed Ssek = ^^y^ at 
any time if the mass of the hole never exceeds Mq. The final entropy of entanglement 
(|6l|) could be made arbitrarily large (as one could see by considering eg. a n-shock 
wave case). 

7. Entropy for 3+1 dimensional black holes 

Let us now turn to the 3+1 dimensional black hole. We will consider only the 
Schwarzschild case. The Bekenstein entropy for this hole is 

5bek = A/A = 47r(2M)V4 = AttM'^ (62) 
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The entropy collected in the form of radiation by an observer at X"*" is also proportional 
to pi, Inl though it requires use of the 'transmission coefficients' T{uj) for its 



computation. 

One is tempted to compare such entropies to the result obtained by carrying out 
the flat space calculation of Srednicki for the case of 3 space dimensions. In the latter 
calculation one considers a massless scalar field, say, in the vacuum state in three space 
dimensions. One traces over the field modes inside an imaginary sphere of radius R, and 
computes the entropy of the corresponding reduced density matrix. In this calculation 
it is convenient to decompose the scalar field into angular modes Yi^m{6,(f)), so that 
we obtain a 1-dimensional problem in the radial co-ordinate r for each such angular 
mode. The different angular modes decouple from each other, so we have to just add 
the entropies resulting from the computation for each mode. The radial co-ordinate 
is taken as a 1-dimensional lattice with lattice spacing a. The entropy is found by 
numerical computation to be [0] 

S ^ .30{R/ay (63) 

Thus we seem to reproduce the ~ dependence expected of the black hole entropy. 
But if we accept that (^) applies to the black hole then we are faced with the question: 
What should be the value of the cutoff a? 

As we now argue, the result (|63D is not the one relevant to the 3+1 dimensional 



black hole. In fact, the 1-dimensional result (|50D is again the relevant one to use. To 
see this, consider decomposing the scalar field in the black hole geometry into modes 
with angular dependence Yi^m{0,(p)- Since we have assumed spherical symmetry, these 
modes decouple from each other. Thus we can consider studying the evolution of 
different fields (labelled by {l,m}) on the 1+1 dimensional geometry obtained by the 
spherically symmetric reduction of the 3+1 dimensional geometry. Proceeding in this 
way, one would need to find the 'entanglement entropy' of fields in one space dimension, 
which (for free fields) is given by the the result (|50|) . 

At this point one sees an important difference between the flat space 3-dimensional 
problem and the 3+1 dimensional black hole. This difference occurs in the number of 
angular modes YJ^m that contribute signiflcantly to the entropy. Consider flrst the flat 
space problem. If we had taken a lattice with lattice spacing a all over 3-space, then 
on the boundary sphere of radius R we could consider angular modes with / ~ R/a. If 
we flrst reduce the Hamiltonian to a sum over modes and then put the r co-ordinate 
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on a lattice, then again it is found that for / >> R/a the degrees of freedom on the 
two sides of the r = R boundary effectively decouple 0. This happens because the 
radial wave equation is dominated at large / by a 'mass term' arising from the angular 
Laplacian, and such a term does not couple neighbouring sites of the r co-ordinate 
lattice. Since —l<m<l, the number of angular modes contributing to the entropy is 
0{{R/aY}, which explains the leading power dependence of the entropy (^) on the 
cutoff a. (Treating the 1-dimensional problem for each angular mode should give a 
ln{R/a) dependence multiplying the dependence ~ {R/aY, but the above argument is 
too crude to distinguish the presence or absence of logarithmic terms.) 

Thus we see that the difference between the a dependence of ( ^OD and (^) can be 
understood in terms of the large number ( 0{{R/a)'^}) of angular modes contributing 
to the entropy in the 3-dimensional flat space problem. But in the 3+1 dimensional 
black hole, we know that most of the radiation comes out in only a few angular modes! 
In fact for a reasonable first estimate of the Hawking radiation one can require that 
only the s-wave modes (/ = 0) emerge from the hole. We now compute the entropy of 
the 3+1 hole with such an approximation. 

While we cannot solve the geometry of the evaporating 3+1 dimensional hole as 
accurately as for the RST model, for the purposes of our calculation we can consider 
the Schwarzschild metric with a time dependent mass M. The surface gravity of the 
black hole is k = (4M)~^. The temperature is T = k,/27i = {87iM)~^. The luminosity 
in the s-wave mode is 

_ dE _ 1 duu _ vrT^ 

' ~ III ~ 2^Jo e-/^ - 1 ~ 12 ^ ^ 

From ( |6^ we compute M{y), the mass of the hole at the T"*" Schwartzschild coordinate 
y. (We take y = at the endpoint of evaporation, so y is negative in the part of X+ 
where the Hawking radiation is being received.) We have 

M(y) = i^r" . (65) 

As mentioned above, we approximate the evaporating geometry by just letting M 
depend on time. Letting v be the Minkowski co-ordinate at T~ . This approximation 
then gives 

dy = -4:M{y)d{ln{vo-v)) (66) 
when the co-ordinate v is close to the value Vq which reflects at r = to move along 
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the event horizon. Integrating (|6^) gives 



ln(fo - Vf) - \n{vo - Vi) = 967r(M; - M/) (67) 

We have Mj = M, and we let Mf be of the order of Planck mass. Further, ln(f j — vq) 
can be ignored compared to ln(fo — Vf). Following the discussion of entropy in the 1- 
dimensional case, we write (for one evaporating field, one 'cut' and rightmovers only): 

S ^ -^ln(a) = -^ln(H (68) 

Here Sv is found from (p6| ) by setting dy = A (for some chosen scale A over which the 
observation of radiation is switched off): 

5v = -^{vo-Vf) ^ (69) 

Substituting in (|6^) we get 

S ^ SttM^ = 25Bek (70) 

We have ignored terms logarithmic in M; these corrections are smaller than the con- 
tribution of the / 7^ modes which we have also neglected. 

We can now compare the result (|70D for the entanglement entropy with the ther- 



modynamic entropy collected at X+. Following [|1^], we compare the change in the 
thermodynamic entropy received at X+ to the change in the Bekenstein entropy of the 
hole: 

_ _ /o°° dxxMx){xey{e^ - 1) - ln(e- - 1)} 
R - dS/dS^^^ - /-rfxx3a(x)/(e^-l) ^^^^ 



Here 



= ^rz,^(^)/[27(a;M)2] (72) 

Lm 



is the absorptivity per unit area of the black hole. In the above used approximation to 
Hawking radiation we used only the I = mode, and let the 'transmission coefficient' 
r be unity for all u. This gives <j{uj) = 1/[27(cijM)^]. Substituting this in ( [TiD gives 
i? = 2, in accordance with (|70D. For more realistic models, taking into a account the 
transmission coefficients Ti^m{^), one obtains i? ~ 1.3 — 1.6 [0. To reproduce the 
effect of nontrivial T{uj) we would need to extend the Srednicki calculation to fields 
with position dependent 'mass term'; we do not investigate this further here. 

We can also compare the above derivations to the direct computation of the entropy 
of the density matrix obtained in the evaporation process; i.e., carry out the calculation 
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analogous to (139|) to (^9|) . We again have S'„ = vr'^T/Se, but now T = (87rM„) ^ = T„ 



Following the evaporation process we find Mn = (n/128e)i/l This gives, as expected, 

"max 7j-3t^ 

^ = E ^ 8^^' = 2^Bek . (73) 

n=l 

In our above application of Srednicki's result, we found that the one space dimen- 
sional formulation was more relevant, rather than the three space dimensional one. 
On the other hand if an observer stands near the horizon of the black hole, she sees 
thermal radiation with power in a large number of angular momentum modes. Then 
it is possible that by carrying out the above calculations with a different choice of 
hypersurface {e.g. with the 'outside' part of the hypersurface corresponding to a static 
frame near the horizon) one would find relevant the equation (|63D. 

8. Discussion 

In this paper we have carried out the computation of Bogoliubov coefficients, stress 
tensor and the entanglement entropy for an evaporating black hole. Before one had 
explicit models of evaporating geometries, such quantities had been calculated only in 
the absence of backreaction. With backreaction, it is possible to obtain for example 
the stress tensor in the region between the event horizon and the apparent horizon. 

Concerning the application of Srednicki's result to the black hole entanglement 
entropy, we have made two assumptions. First, we have assumed that the infrared 
divergence comes from modes with wavelength very large compared to the system 
size; after these modes are removed, the entropy can be split into a contribution from 
rightmoving modes and a contribution from leftmoving modes. Second, we assumed 
that when we cut the region R out of a line, the entropy of system, say, can be split 
into two parts, one coming from each of the two cuts at the two ends of R. What we 
do now is offer some heuristic arguments to justify these assumptions. 

First we wish to understand the appearance of the logarithmic dependence in the 
entanglement entropy. Consider a segment of the real line, < x < Ii + I2, divided 
into two regions near x = Ji by a 'cut' of length a << Ii, I2. Further, assume Ji < I2. 
The scalar field we take to vanish at x = 0, x = Ji + 12- What is the entanglement 
entropy of Ji with I2 (with cutoff scale a)? 

Suppose Ii = 12- The field modes have wave numbers 

/c = — — , n = l,2,... . (74) 

+ ^2 
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The mode = 7r/(/i + 12) is 'shared' between the two sides /i,/2, and we assume that 
it gives an entropy s. Now consider modes with 1 ~ n ~ (Ji + l2)/a. For any scale oj~^ 
of the wavelength we can make localized wavepackets just as was done in section 6 (eq. 
(p^)). The number of such wavepackets is ~ n. Only the wavepacket that overlaps 
both sides of the cut at x = Ji contributes to the entanglement entropy; and we again 
take this entropy to be s. (This is an assumption.) Equivalently, we find that each 
mode (|7^) contributes ^ to the entropy, which becomes 

S^sY-^sln-. (75) 
tin a 

Now suppose instead that Ji << I2. For k^^ » Ii, the mode essentially vanishes 
over Ji, and so cannot entangle this region with I2. For wavelengths ^ Ii we make 
localized wavepackets in the region < x < ax, a ~ 1, just as in the discussion above, 
and thus find again eq. (^), where we note that Ji is the smaller segment. 
Now we address a more complicated case. We have the segments 

• Ji : < a; < Ji 

• R: h < X < I1 + R 

• h: I1 + R< X < I1 + I2 + R . 

Let R « Ii « I2. h and R are separated by a cut of size ai << R, and R 
and I2 are separated by 02 << R. We want the entanglement entropy of R with the 
remainder {i.e., Ii U I2). Again, we do not need to consider modes with wavelength 
» Ii, since these effectively vanish over R, Ii. The modes with R < k~^ < Ii give 
a contribution S" ~ sln^. The modes with k~^ ~ R can be made into wavepackets 
that get 'partitioned' at two places; x = Ii,x = Ii + R. These give the contributions 
S^sln—, S^slia—. Overall, we then get 

ai ' a2 ' ° 

R R I 

S ^ s\n hsln hsln-i . (76) 

ai a2 R 

If oi = 02 = a, Ii = I2 = I, we can write 

5~2sln- + sln- , (77) 
a R 

which resembles (|50D with k,2 = The value of s we can fix by the direct calculation 

of the entropy in the 1+1 dimensional black hole evaporation preceeding eq. (ff3). For 
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the modes with ~ i?, one can clearly make a breakup between right and left movers. 
Thus we conclude s = |, after doubling up the obtained answer for the rightmovers 
alone. 

The above gives a heuristic understanding of the entropy of entanglement, which 
should be useful in applying the result with a variety of boundary conditions. 

In conclusion, we note that the 'exponential expansion' of coordinates near the 
horizon gives rise to thermal radiation, as was shown by Hawking. By using the result 
of Srednicki, the same coordinate transformation gives the entanglement entropy pro- 
duced by the thermal radiation. Thus we seem to be one step closer to understanding 
the nature of black hole thermodynamics. 
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Figure caption 

Figure 1. The black hole geometry formed by an incoming shock wave (thick 
line with arrows). Evaporation occurs in region II; regions I and III are linear dilaton 
vacuua. 
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